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Abstract 



We determine the chiral condensate from simulations of quenched lattice QCD with 
Wilson fermions. Our measurements have been obtained with high statistics at three 
values of the gauge coupling, corresponding to UV cutoffs in the range 2 — 4 GeV. 
Several improvements have been made with respect to earlier lattice computations. 
The most important are the non-perturbative renormalization of the condensate, the 
use of the tree- level improved Clover action and the reduction of the systematic error 
due to uncertainties in the lattice calibration. Our result for the chiral condensate 
in the MS scheme is 

(^) MS (Ai= 2 GeV) = - 0.0147(8)(16)(12) GeV 3 = - [245(4)(9)(7) MeV] 3 
where the first error is statistical, the second is due to the non-perturbative renor- 
malization and the third due to the lattice calibration. 
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1 Introduction 



QCD has a remarkably small mass gap. Its pattern is described by a nearly 
degenerate isospin multiplet containing pseudoscalar mesons, but not their 
partners of opposite parity (vector mesons). This implies that chiral symmetry 
is spontaneously broken in massless QCD. This should be signalled by the 
non-vanishing of some order parameter (for a given order parameter this is 
a sufficient but not necessary condition for symmetry breaking). The chiral 
condensate is such an order parameter. The standard expectation is that it 
does not vanish in the chiral limit. However, it cannot be a priori excluded (see 
ref. [1] and, more recently, [2]) that it might tend to zero in this limit, while 
chiral symmetry breaking is nevertheless realized through the non- vanishing of 
the pion decay constant. A more general option is that the chiral condensate 
has a small non-zero value [2]. 

The intrinsically non-perturbative nature of the chiral condensate implies 
that its determination from first principles is not straightforward. QCD sum 
rules and lattice computations have been used in order to estimate its value. 
The lattice regularization is particularly well suited for the computation of 
low-energy non-perturbative quantities from first principles. Although lattice 
measurements are affected by a series of approximations (finite configuration 
ensemble, finite volume and lattice spacing, non-zero quark mass, quenched 
approximation), they can be controlled and systematically improved, at least 
in principle. Two of these sources of error have the most negative impact on 
the credibility of the lattice determination of the chiral condensate, namely 
the quenched approximation and the chiral extrapolation of data computed at 
non-zero quark masses. Quenching is destined to stay with lattice computa- 
tions for some time, since realistic simulations with light sea quark masses are 
still unaccessible to present-day computers. Quenched simulations with fairly 
light valence quarks are more accessible, but the zero quark limit is elusive for 
hadronic quantities^. 

Interestingly, the two approximations (quenching and chiral extrapolations) 
are interwoven by the presence of quenched chiral logarithms in many phys- 
ical quantities, such as the chiral condensate and the pion mass (but not in 
the pion decay constant, in the isospin limit) [6]. There has been some debate 
recently on the evidence for quenched chiral logarithms in lattice results ob- 
tained with light quark masses [7]. However, for the bulk of quenched lattice 
results, obtained at the strange quark mass region, the power dependence on 

1 Note however, that it is possible to measure the chiral condensate with staggered 
fermions at zero quark mass, using the Banks-Casher formula [3] (see for example 
refs. [4]). More recently, the Shrodinger Functional formalism has been used for 
Wilson fermions in order to obtain lattice renormalization constants at zero quark 
mass [5]. 
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the quark mass dominates the logarithms. In such simulations many physical 
quantities known from experiment, including those related to the chiral con- 
densate (e.g meson masses, decay constants), have been determined with a 
5-15% accuracy (see for example ref. [8]). Also, the recent lattice determina- 
tions of the quark masses [9] agree with the sum rules predictions. The good 
agreement of many quenched results with experiment 0, obtained in a quark 
mass region where chiral singularities can be ignored, implies that the regular 
dependence of these observables on the quark mass is adequately reproduced 
in the quenched approximation. 

Given the above observations, we have performed an extensive quenched lattice 
computation of the chiral condensate, with Wilson fermions, using the data 
accumulated by the APE collaboration in the last few years. The valence quark 
masses used in these simulations are in the range of the strange flavour. Several 
important improvements have been possible compared to earlier computations: 

• We have improved our statistics by computing the condensate on ensem- 
bles of 0(100) gauge field configurations. Most previous computations with 
Wilson fermions were performed on O(10) configurations. 

• Reliable measurements have been performed at two values of the lattice 
gauge coupling g%, corresponding to UV cutoffs of ~ 2 GeV and ~ 3 GeV, 
using both the Wilson and the tree-level SW-Clover[10,ll] actions. Results 
obtained with the Wilson action at fixed UV cutoff (i.e. fixed lattice spacing 
a) suffer from 0(a) discretization errors; those obtained in the tree- level 
Clover formalism suffer from 0(agQ) errors. Thus we can study the influence 
of finite lattice spacing effects. We also present results at a third value of 
the gauge coupling (corresponding to an UV cutoff of ~ 4 GeV), which are 
however somewhat unreliable, due to the small physical size of the lattice. 

• Different determinations of the chiral condensate (expressed in lattice units) 
have been implemented. They all derive from the same Ward Identity (WI), 
obtained from the lattice action, and are equivalent only in the chiral and 
continuum limits. They are implemented, however, at non-zero lattice spac- 
ing and quark masses and subsequently extrapolated to the chiral limit. 
Thus, by comparing the chiral condensate from different determinations, 
we can cross-check the reliability of the extrapolations. 

• The systematic errors arising from the chiral extrapolation and the un- 
certainties in the determination of the lattice spacing are amplified in the 
standard computation of the chiral condensate. The chiral extrapolation am- 
plifies the error because we extrapolate from data obtained in the strange 
quark region. The uncertainty in the determination of the lattice spacing is 
amplified because the condensate is a dimension-3 matrix element. We have 
obtained the lattice condensate in physical units in a way that avoids these 



2 There are several exceptions, such as the disagreement of the experimental and 
the lattice result (in the quenched approximation) of the decay constant ratio fx/ fn- 
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problems. 

• The chiral condensate measured on the lattice at a given fixed lattice spac- 
ing is a bare quantity which must be renormalized. In previous studies, 
its renormalization constant was calculated in lattice Perturbation Theory 
(PT) at 1-loop. We will show that the perturbatively renormalized chiral 
condensate suffers from large 0(qq) corrections, which render its determi- 
nation unreliable (an identical conclusion was drawn for the quark mass in 
ref. [9]). In the present work, we have used non-perturbative (NP) estimates 
of the required renormalization constant, which are known in the Regular- 
ization Independent (RI) renormalization scheme (also known as the MOM 
scheme). This removes lattice PT completely from the renormalization pro- 
cedure. We only use continuum PT in order to transform the results from 
the RI to the intrinsically perturbative MS scheme. Also, RG running has 
been done in PT in order to express all results in the conventional scale of 
2 GeV. All PT coefficients used in this work are at NNLO. 

Our best estimate for the chiral condensate agrees with the earlier lattice 
determinations of refs. [12,13] with greatly improved statistical accuracy. For 
the first time we are also able to estimate the systematic uncertainties in a 
reliable way. Our result is also in excellent agreement with the one obtained 
from sum rules (see [14] and references therein). 

We now outline the contents of this paper. In sect. 2 we fix our notation and 
give the basic definitions. In sect. 3 we discuss issues related to our renor- 
malization strategy. In particular we review the most important properties 
of the NP renormalization of the lattice condensate (see Appendix A for de- 
tails). We also review the connection between the RI and MS renormalization 
schemes (details can be found in Appendix B). In sect. 4 we discuss how 
the chiral condensate can be determined from the lattice regularization with 
Wilson fermions. Chiral symmetry, explicitly broken by the Wilson fermionic 
action, is recovered in the continuum and chiral limits through the implemen- 
tation of lattice Wis; see ref. [15]. With the aid of such a WI we obtain the 
proper lattice definition of the chiral condensate which, for Wilson fermions, 
is subject to a power subtraction. Moreover, from this WI we derive the well 
known Gell-Mann-Oakes-Renner (GMOR) relation [16] for the condensate. 
Some technicalities related to the Wis with Clover fermions are discussed in 
Appendix C. In sect. 5 we present our results and compare them to those of 
previous computations of the chiral condensate with Wilson fermionsQ Our 
conclusions are summarized in sect. 6. 



3 No comparisons were made with results obtained with Kogut-Susskind fermions, 
which are usually presented unrenormalized. We also point out that this work is 
limited to zero-temperature QCD. For a recent review on lattice studies of chiral 
symmetry restoration at finite-temperature QCD see ref. [17]. 
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2 Definitions and generalities 



In this section we define the quantities of interest to this paper, in order to 
fix our notation. We work in the lattice regularization scheme proposed by 
Wilson [18]. The partition function of the theory is 

Z = f VipvjVUn exp (S g - S f ) , (1) 



where the gluonic action is given by 

6 



Sg ~ — 

9o P 



£{l-^r[C/p + [/t,]} (2) 



and the fermionic action by 



x,/i la 

+a 4 E^(^)(M + -)V(x) . (3) 

x a 

In standard notation, a is the lattice spacing, ip{x) is the quark field (flavour 
indices are implied), Up is the Wilson plaquette, U^(x) is the lattice gauge 
link and g Q is the bare coupling constant. The diagonal bare mass matrix is 
denoted by M . In most of this work we will assume mass degeneracy; i.e. 
all elements of M are equal to m . As is well known, for Wilson fermions, 
the quark mass, besides a multiplicative renormalization, is also subject to an 
additive one. Defining m = m — mc, the renormalized is 

m R = Z m m = Z m [m - m c ] (4) 



with Z m the multiplicative renormalization constant. The chiral limit is then 
m — > 0; i.e. m — > mc- At tree-level, mc = — 4/a. 

We define the scalar and pseudoscalar densities as: 



S'(x)=$(x)^(x) 
f - ^ 

P J (x) =i)(x)^—i)(x) (5) 
and the vector and axial currents as 
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f - xf 

A l{x) = i)(x)^^ b —i)(x) , (6) 

where are the generators of the SU(Nj) flavour group in the fundamental 
representation (e.g. Pauli matrices for Nf = 2, Gell-Mann matrices for Nf = 3 
with / = !,..., N* — 1). These generators satisfy: 




(7) 



where ff gh are the SU(Nf) structure constants, d^ gh are symmetric coefficients 
and A represents the Nf x Nf identity matrix. We will also consider the 
corresponding singlet bilinear scalar operator 



S°(x) = %fj{x)\ ip{x) = ip(x)tp(x) 



The operators defined above are to be understood as bare lattice quantities^. 
They are subject to multiplicative renormahzation. The renormahzation of a 
generic bilinear operator 0[ = tjjT(X^ /2)ijj is given by (r stands for any Dirac 
matrix) 

6r(g%, m R , n) = lim[Z (^o, ^)0 T {gl, m, a)] , (9) 



where gu is the renormalized gauge coupling and \i the renormahzation scale. 
Note that in the bare operator the bare mass mo has been traded off for the 
more convenient (from the point of view of the chiral limit) subtracted mass 
m. Renormalized quantities will be denoted in general by "hats"; e.g. the 
renormalized quark propagator is 

S{p) = ^n[Z q S(p)} , (10) 



where Z^ 2 is the quark field renormahzation. 

4 This means that all integrals J d 4 x which appear below, are to be interpreted as 
sums (a 4 J2 x ) which run over all lattice sites, labelled by x; 5(x) stands for a 5 x q, 
etc. The use of integrals instead of sums and Dirac functions instead of Kronecker 
symbols is a question of notational convenience and should be transparent. 
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The properties of the renormalization constants of the operators of eqs. (5) and 
(6), obtained in refs. [15] (for a recent treatment see [19]), are as follows: The 
renormalization constants Z$ and Zp of the scalar and pseudoscalar densities 
are logarithmically divergent; thus they depend on the coupling g% and the 
renormalization scale in lattice units a\x. Those of the vector and axial currents 
[Zy and Za respectively) are finite; thus they depend only on the coupling g\. 
A mass independent renormalization scheme is implied. 



3 Renormalization scheme dependence of the chiral condensate 



In this section we discuss some issues related to the choice of renormalization 
scheme for the chiral condensate. The renormalization procedure which gives a 
finite condensate starting from its bare (lattice) counterpart is postponed until 
the next section. Here we will concentrate on how to connect the condensate 
in the RI renormalization scheme, obtained non-perturbatively, to the one in 
the MS scheme. This is necessary in view of the fact that the RI scheme is 
most convenient for lattice computations, as discussed below, whereas the MS 
scheme is conventionally preferred when expressing "final" physical results. 

Let us sketch step by step the renormalization procedure for the chiral con- 
densate. The breaking of chiral symmetry by the Wilson term of the fermionic 
action Sf induces some subtleties in the renormalization of the lattice conden- 
sate. These have been dealt with in ref. [15] and will be fully exposed in sec. 4. 
Here we simply anticipate the basic result: by using lattice Wis and by re- 
quiring that in the continuum and chiral limit the nominal Wis are recovered, 
we obtain that the correct renormalization of the lattice chiral condensate is 
given by 

{W){H) = Z P (fia) \{W){a) - 6ol , (11) 



where (ipip)(a) = (S°) stands for the bare chiral condensate, computed non- 
perturbatively on the lattice, and b is a cubicly diverging subtraction. Note 
that Zp is the renormalization of the non-singlet pseudoscalar density. Its de- 
termination from 1-loop lattice perturbation theory (see refs. [20-22] for the 
perturbative renormalization of quark bilinears) suffers from large uncertain- 
ties due to the presence of "tadpole" diagrams [23]. In order to avoid this 
problem non-perturbative renormalization techniques have been developed 
[5,12,24,25,40] consisting in computing the renormalization constants at fixed 
UV cutoff (ranging, in present-day simulations, in the region a" 1 ~ 2—4 GeV). 
This is a tradeoff between higher order corrections (present in PT) and lattice 
artifacts (present in all non-perturbative methods). The latter, however, can be 
systematically controlled by improving the action and operators in the spirit of 
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Symanzik's work [5,10,11,26]. Thus in principle non-perturbative evaluations 
of the renormalization constants are preferable. The finite renormalization 
constants (Z v , Z A and the ratio Z P /Z S ) can be obtained non-perturbatively 
with the help of Wis (see [5,12,15,19,24,40]). All renormalization constants 
(including the diverging ones Zs and Zp) can be obtained non-perturbatively 
with the help of the more recent "Non-Perturbative" (NP) method of ref. [25]. 
Results of the NP renormalization constants of quark bilinear operators can 
be found in refs. [25,27,28]. 

In practice, the necessity of non-perturbative renormalization was first demon- 
strated in the measurement of several B-parameters of the AS = 2 four- 
fermion operators; see refs. [29]. More recently, it has also been shown in 
ref. [9] that the discrepancy, observed in the past between quark masses ex- 
tracted from the Vector and Axial Wis, was due to the very poor determination 
of Z P in PT. Use of the NP value led to an excellent agreement between the 
two determinations. The same conclusions are drawn in the present work for 
the chiral condensate. 

We will therefore use the NP value of Zp in the present work. The underlying 
theory and method of computation has been fully exposed in [25] ; it consists in 
applying the RI (or MOM) scheme on the lattice. The amputated Green func- 
tion of the insertion of the operator P in the quark propagator is computed at 
fixed external momentum in the deep Euclidean region. The Green function is 
subsequently projected by a "suitable projector" and the RI renormalization 
condition fixes the renormalized projected Green function be equal to its tree- 
level value at a given scale \x. This condition can then be solved for Zp; the 
result is a fully non-perturbative estimate of the renormalization constant at 
fixed UV cutoff. For this procedure to be reliable \x must satisfy the conditions 
/i <C Oia) (in order to avoid discretization errors) and /i ^> Aq C d- The former 
bound is approximately satisfied in practical simulations; we will be using the 
renormalization constants of ref. [28] at the scale \i ~ a -1 . The latter bound 
is necessary for two reasons, (i) The RI scheme imposes that the renormalized 
operator transform as an irreducible representation of the chiral group. This 
is true only at large scales, where chiral symmetry breaking effects are negli- 
gible, (ii) At large scales we avoid higher order corrections in the continuum 
perturbative expansion in which the RI-MS matching coefficient is calculated; 
see below. Unlike MS, the RI renormalization condition does not depend on 
the regularization chosen, hence its name, Regularization Independent (RI). 
The explicit formulae of the RI renormalization condition and a discussion 
of some subtleties related to wavefunction renormalization can be found in 
Appendix A. 

Thus the chiral condensate of eq. (11), renormalized with the NP method 
outlined above at scale ji ~ a -1 , is expressed in the RI scheme. In order to 
enable comparisons with other methods, we must convert our result to the 
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MS scheme^. This is expressed in terms of a perturbative finite matching 
coefficient 

(W>) m (j*) = kZ m > m {^) m {ii) . (12) 



The matching coefficient AZ MS ' RI , known in PT to NNLO (see refs. [28] 
and [30]), is given in Appendix B. Note that the renormalization constants 
obtained in the RI scheme are in general gauge dependent (we work in the 
Landau gauge). This dependence is cancelled by the gauge dependence of 

In order to express our result at the conventional scale of p! = 2 GeV, we use 
perturbative RG running 

W) M V) = y^(^) MS (/i) • (13) 



Finally, it is convenient to express the chiral condensate in a scale independent 
way. We define the RGI chiral condensate as 

<^> RGI = - Jrr • (14) 
cf s (/i) 



This quantity does not depend neither on the scale fi nor on the renormal- 
ization scheme. The evolution coefficient c]y^(//) is known to NNLO; see Ap- 
pendix B. 



4 The chiral condensate from lattice Wis with Wilson fermions 



In this section we obtain the proper definition of the chiral condensate with 
Wilson fermions from lattice Wis. The subject of chiral symmetry breaking 
by the Wilson term and its restoration in the continuum limit with the aid of 
Wis is well understood (see the original works [15], a recent review [19] and 
a recent theoretical treatment of the subject [31]). Thus we only present in 
some detail those results essential to the understanding of the lattice deter- 
mination of the chiral condensate. First we shall briefly digress to vector Wis 
in order to demonstrate that the quark mass renormalization is the inverse 
of the renormalization of the scalar density. Then we will obtain the proper 

5 The continuum perturbative calculations related to the MS results quoted 
throughout this work have been performed in the NDR regularization scheme. 
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definition of the lattice chiral condensate and its renormalization using axial 
Wis from which the Gell-Mann-Oakes-Renner relation [16] will be derived on 
the lattice. 



The local SU L (Nf) x SU R (Nf) chiral transformations of the fermionic fields 
are defined as follows: 



Si/j(x) = i 
5?p(x) = —iijj(x) 



ft n A/ ft ^ 



ip(x) 



' ff ^ fr \ Xf 



(15) 



With degenerate bare quark masses tuq, the global vector transformations 
are a symmetry of the action. From the corresponding local transformations 
(eqs. (15) with a^ A = 0) vector Wis can be derived, which indicate that the con- 
served lattice vector current is a point-split operator. The no-renormalization 
theorem applies to this operator. On the other hand, the lattice vector local 
current of eq. (6) is subject of a finite renormalization Zy^gfy ^ 1, due to 
its non-conservation on the lattice. A detailed treatment of these results can 
be found in refs. [15,19]. Here we illustrate the point by writing down a typ- 
ical lattice vector WI. For simplicity, the flavour / of the vector variation of 
eq. (15) is chosen to be non-singlet and non-diagonal (e.g. the lowering oper- 
ator \f = A 1 — i\ 2 ). The two quark masses related to this flavour are denoted 
by mi and m 2 (i.e. they are non-degenerate). Then we obtain the WI 



Sa v (x 



•<K?(o)> = o 



Z v J2W x {Vf(x)V°(0)) = (m 2 - mi )(Sf(x)V9(0)) , 



(16) 



where aV%f(x) = f(x) — f(x — fi) is an asymmetric lattice derivative. We 
have imposed x ^ 0, in order to avoid contact terms. We can now multiply 
both sides of the WI by Zy and require that the renormalized quantities 
obey the nominal vector WI. This implies that the product of the quark mass 
and the scalar density is renormalization group invariant. In other words, the 
renormalization of the scalar density S(x) is the inverse of the multiplicative 
renormalization of the quark mass: 



(17) 



Far less straightforward is the implementation of the axial symmetry with Wil- 
son fermions, because of the presence of the chiral symmetry breaking Wilson 
term in the action. The basic idea is that, by imposing suitable renormal- 
ization conditions, PCAC is recovered in the continuum. Consider the axial 
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WI, obtained from eq. (15) with ay = 0, arising from the variation of the 
non-singlet pseudoscalar density. For /, g ^ we have: 



5a A (x) ^ ^ 



1 

-(^(x)P"(0)> 



— <S°(0)> = J2^(Ai(x)P 9 (0)} - 2m (pf(x)P*(0)} (18) 



The l.h.s. of this equation is the variation of the operator which vanishes when 
x 7^ 0. The r.h.s. is the variation of the action from which the lattice version 
of the standard PCAC relation is recovered. The operator X* in the above 
WI arises from the variation of the Wilson term under the axial transforma- 
tion 6 . Unlike the vector current case, X- f cannot be cast in the form of a 
four-divergence. It is a dimension-4 operator which, in the naive continuum 
limit vanishes, being of the form (axdimension-5 operator). However, it has 
divergent matrix elements beyond tree-level. Its mixing with operators of equal 
and lower dimensions can be expressed as follows [15]: 

X f (x) = X f (x) + 2mP f (x) + (Z A - l)V£Aj(a;) , (19) 



where naive dimensional arguments tell us that the mixing constant Z A {g^) is 
finite, whereas Tn(gQ,mo) diverges linearly like a™ 1 . Logarithmic divergences 
have been shown to be absent at all orders in PT [32]. The two mixing con- 
stants Z A and m, and consequently X (x), are defined by requiring that the 
renormalized axial current satisfies the nominal (continuum) axial WI. This 
requirement is satisfied upon demanding that, in the continuum limit and for 
vanishing quark mass, the correlation functions of X 9 vanishes except for lo- 
calized contact terms. For the insertion of X in correlation functions such as 
(X (x)P 9 (0)) the contact terms, by flavour symmetry, must have the form 



(X f (x)P 9 (0)) = -^-b 5(x)5 f9 + hS^USix) , (20) 
N f 

where the last term on the r.h.s. is a Schwinger term (it vanishes upon inte- 
gration over / d A x). Combining eqs. (18), (19) and (20) we arrive at 



Z A V»(Al(x)P 9 (0)) -2(m -m){P f (x)P 9 (0)) - M /9 D<5(a;) (21) 

6 The original WI, obtained as described above, involves a point-split axial current. 
This current, expanded in powers of the lattice spacing a, is equal to the local 
current of the above WI plus terms which are absorbed in a redefinition of (see 
refs. [15,19] for details). 
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= -5(x)5f±-[(S°)-b 

By multiplying both sides of the above by Zp, we obtain a completely finite 
expression, which is the renormalized WI in the continuum. The &i-subtraction 
on the l.h.s. compensates the contact terms arising in the correlation functions 
V^(Aj l (x)P 9 (0)) and (pf \x)P g (0)) , when x is close to the origin. Otherwise, 
we have recovered the nominal WI, which is now satisfied by renormalized 
operators. This implies the following useful results: 

(1) The product = Z^A? must be interpreted as the renormalized axial 
current (Za(<7q) 7^ 1 is a finite renormalization). 

(2) The chiral limit is to be defined as the value mc of mo for which the 
difference mc — m(gQ,m c ) vanishes. We see from eq. (21) that the product 
(mo — m)Pf is renormalization group invariant. Thus, we can write the renor- 
malized quark mass as: 

m R = Z m [m - m(m )] , (22) 

and obtain for the renormalization of the pseudoscalar density 

Z P = Zl] . (23) 



Combining eqs. (4), (17), (22) and (23), we obtain for the ratio Z P /Z S 

Z P m — m(m ) 
Z s m - m c 



(24) 



(3) Two useful expressions are derived from eq. (21), when 1 / (i.e. all 
terms proportional to S(x) vanish). The first expression consists simply in the 
definition of the ratio (for / = g): 

_ 2(m -m) fd*xVl(4(x)Pf(0)) 
' Z A " Jd 3 x(Pf(x)Pf(0)) ' 1 ' 



where the spatial derivatives of the divergence of the axial current vanish 
under the integration. This quantity, once Za is known, can be used for the 
determination of (mo — m). The second useful expression arises upon applying 
the LSZ reduction formula: 

Z A V%{0\Al(0)\P) = 2(m -m)(0|P / (0)|P) , (26) 



where \P) is the lowest lying pseudoscalar state (in the chiral limit this is the 
Goldostone boson). Both equations are lattice versions of the standard PCAC 
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relation. 

(4) Last but not least, the ^-subtraction on the r.h.s. of eq. (21) implies that 
the proper definition of the lattice chiral condensate with Wilson fermions, 
determined by axial Wis, is given by 



(W)«*=\(S°)-bo 



(27) 



The above "subtracted chiral condensate" is a logarithmically divergent quan- 
tity, which is multiplicatively renormalized by Z P : 



<H>) 



(S°) - b 



(28) 



To obtain an expression from which the chiral condensate can be computed in 
a convenient way, we need to integrate both sides of WI (21) over all space- 
time. Upon integration the Schwinger ^-term on the l.h.s. vanishes. For the 
remaining terms on the l.h.s. two possibilities arise: either work directly in the 
chiral limit or work with a non-zero quark mass and subsequently extrapolate 
to the chiral limit. If we work directly in the chiral limit, the (m — m) term 
vanishes. The integrated V^A£ would also vanish, being the integral of a 
four-divergence, if chiral symmetry were not broken (absence of Goldstone 
bosons). In this case the WI implies a vanishing chiral condensate. If the 
symmetry is broken, the presence of a Goldstone boson guarantees a non- 
vanishing surface term upon integrating V^A^. Thus, also the r.h.s. of the 
WI (i.e. the chiral condensate) is non-zero. In the present work, in accordance 
with standard numerical computations, we will be working with a finite quark 
mass subsequently extrapolated to the chiral limit. Then the surface term of 
the divergence of the axial current vanishes, and the integrated WI becomes 
(for f = g) 

-L(^) sub = hm 2(m - m) [ d*x(P f (x)P f (0)) . (29) 

Nf mo^mc J 



From the above expression, we can also derive two other useful identities for the 
chiral condensate. They are equivalent, in theory, to the one above, in the chiral 
limit and up to discretization errors. Since in practice they are implemented at 
finite quark mass and lattice spacing, they may reveal the importance of these 
systematic effects. Both expressions are obtained by writing the correlation 
function [P* (x)P* (0)) as a time-ordered product and by inserting a complete 
set of states in standard fashion. The spatial integration / d 3 x projects zero- 
momentum states. Contributions from higher mass states vanish in the chiral 
limit. Upon performing the time integration / dt, we find 

-L(^) sub = - hm (m °^ M) (0|P(0)|P) 2 , (30) 
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where rap is the mass of the pseudoscalar state \P). Note that the pseudoscalar 
operator without a flavour index is defined above as, say, P = u^^d (and the 
pion state has a du flavour content). The same is true for the axial current 
(without flavour index) = U'j^d. 

Another expression for the chiral condensate is obtained by substituting in the 
above equation the matrix element (0\P(0)\P) by Vg (0|A M (0)|P) (see eq.(26)). 
The matrix element of the axial current is parameterized in terms of the pion 
decay constant: 

(0\A^x)\P(£)) = ifpp^eM-ipx) , (31) 

(fp is the decay constant of the pseudoscalar state \P)) and thus 

Vft<0|A M (0)|P(0)) = /pm!. (32) 

Combining eqs. (26), (30) and (32) we find 

T7-(^)«i6 = - lim A / Pmp —, , (33) 

which is the familiar Gell-Mann-Oakes-Renner (GMOR) relation [16]. Note 
that the non-vanishing of the chiral condensate in the above equation implies 
the well known linear dependence of the pseudoscalar mass squared from the 
quark mass. 

All determinations of the chiral condensate discussed above are subject to 
taking the chiral limit (m — > m c ) and the continuum limit (a — > 0). In 
practice we work at non-zero quark mass ttiq (fixed hopping parameter k) and 
lattice spacing a (fixed inverse square coupling f3). At fixed (3 the chiral limit 
is normally taken by extrapolating the lattice data to kq ( see comments in the 
introduction on the reliability of this extrapolation and the presence of chiral 
logarithms). Assuming that this extrapolation is reliable, we obtain a result 
which is subject to systematic errors primarily due to the finiteness of the 
lattice spacing. These effects are 0(a) for Wilson fermions. To reduce these 
errors, improved actions in the spirit of Symanzik's programme [26] must be 
used. We have used the tree-level improved Clover action and operators of 
refs. [10,11], which leaves us with 0(ag^) discretization errors. Comparison 
of results obtained with the Wilson and Clover actions enable us to estimate 
finite lattice spacing effects. Some technical details on the implementation of 
chiral Wis with the tree- level Clover action are given in Appendix C. 
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5 Lattice measurements of the chiral condensate 

Having shown how lattice Wis uniquely determine the subtracted chiral con- 
densate with Wilson fermions, in this section we will explore several ways of 
measuring it from simulations. We will first discuss some lattice details in sub- 
sect. 5.1. In subsect. 5.2 we will explore several standard methods, previously 
used in the literature, for extracting the chiral condensate from lattice simu- 
lations. They are based on eqs. (29), (30) and (33). Their main shortcoming 
is that, upon passing from lattice to physical units, the systematic error of 
the condensate due to the uncertainty in the lattice calibration is amplified. 
Another (minor) error amplification is due to the extrapolation of the result 
from the strange quark region to the chiral limit. In subsect. 5.3 we will present 
what we consider the best method of measuring the lattice condensate, which 
sidesteps these problems. The reader uninterested in lattice technicalities may 
skip these subsections. Finally, in subsect. 5.4 we discuss our final results. 

5.1 Computational details 

In this work we have used the results of the best runs performed by the APE 
group in the last few years. Both the Wilson (W) and the tree level Clover 
(C) actions have been used at three values of the coupling (3 = and 
several quark masses in the strange quark region (the range of the lowest-lying 
pseudoscalar meson in these runs is in the range — 2mjc). Table 1 shows 
the parameters of the runs from which we have extracted the necessary masses 
and matrix elements. A discussion of the technical details can be found in [8] 
and [9]. Here we just state that all statistical errors have been estimated with 
the jacknife method by blocking the data in groups of 10 configurations. We 
have checked that our error estimates are roughly independent of the blocking 
size. The APE collaboration has also performed dedicated runs at the same 
(3 values with both the Wilson and the tree-level Clover actions, in order to 
obtain NP estimates of the renormalization constants in the RI scheme [28]. 
The parameters of these runs and the results for the renormalization constants 
of interest to the present work are reported in table 2. More results and details 
on the error quoted can be found in ref . [28] . The main qualitative conclusion 
of ref. [28] (as seen in Table 2) is that the NP method gives excellent estimates 
of the renormalization constants Za and Z$, whereas Zp (and subsequently 
the ratio Z P /Z S ) are subject to a bigger systematic error. 

In lattice simulations we use dimensionless quantities; we thus define the op- 
erators (in lattice units)Q 

7 The lack of an 1/Nf factor in (xx) simply reflects the fact that, since we work in 
the quenched approximation and with degenerate quark masses, there is no flavour 
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Table 1 

Summary of the parameters of the runs with the Clover (C) and Wilson (W) fermion 
actions used in this work. We also show the time intervals in which correlation 
functions were fitted in order to extract meson masses and decay constants. The 
bottom line is our preferred value for the lattice spacing, obtained by fixing the 
mass of the vector X*-meson. The error is statistical. 

V{x) = a 3 P(x) 
A II (x) = a 3 A ll (x) 



{xx) = a z —{^)sub ■ 



(34) 



The mass and decay constant of the pseudoscalar meson in lattice units are 
denoted by Mp = amp and Fp = afp. The extraction of quantities of interest 
in lattice units at fixed quark mass (e.g. Mp, Fp, ap) consists in fitting the 
large time behaviour of correlation functions of suitable operators. Since this 
is a standard procedure we will not present it here; the interested reader can 
find all relevant details of our analysis in ref. [8]. We recall that in terms of 
the Wilson hopping parameter k, the bare quark mass is given byQ 



1 1 

a(m - m c ) = - h 0(a) . 

2k Zkc 



(35) 



To pass from lattice units to physical ones, the experimental value of several 



distinction in our numerical simulations 

8 Note that eq. (35) is also valid in the Clover case; i.e. even in the Clover formu- 
lation it suffers from 0{a) corrections. The corresponding C(a)-improved relation 
will not be implemented as it subject to some complications of a technical nature, 
which we fully explain in Appendix C. 
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83(21 
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85(2") 


0.72(1) 


85(2") 


74(11 


Zp 


0.41(6) 


0.45(6) 


0.47(5) 


0.50(5) 


0.55(3) 


0.57(4) 


ZpjZ s 


0.49(6) 


0.66(8) 


0.56(5) 


0.69(7) 


0.65(2) 


0.77(5) 


Z A 


1.05(3) 


0.81(1) 


1.02(2) 


0.81(1) 


1.01(1) 


0.82(1) 



Table 2 

Parameters of the runs used for the NP calculation of the renormalization constants 
and values of Z$, Zp, Zp/Zs and Za at a scale fxa ~ 1. All results are in the chiral 
limit. The error includes both statistical and systematic effects. 

physical quantities is required (one for setting the scale and one per quark 
flavour for setting the corresponding quark mass). In practice more than one 
physical quantity is used for the determination of the lattice spacing, in order 
to estimate the systematic error of the lattice calibration. In the present work 
we have set the scale a -1 from the vector meson mass Mx* = 0.8931 GeV. As 
suggested in ref. [8], this gives the most reliable estimate, since it corresponds 
to the strange quark mass region, where our simulations have been carried 
out, and it is free of any uncertainties due to renormalization. The result (ob- 
tained with the methods of ref. [8] on the dataset of ref. [9]) is shown in Table 
1. We see that the statistical error is about 2%-5%. The systematic error on 
the lattice spacing has been estimated as the spread of the a -1 values obtained 
from 10 other physical quantities^], as listed in Table 11 of ref. [8]. This er- 
ror is about 7%. Thus, we quote an overall 8% error due to the uncertainties 
in the lattice calibration. The discrepancies between the various lattice cal- 
ibrations is believed to reflect, at least in part, the effect of quenching. We 
will also be using some physical quantities in order to determine the chiral 
limit, quark masses etc. In the region of the up and down quark masses we 
have used the experimental values = 0.138 GeV and /„- = 0.1307 GeV. 
In the region of the strange quark mass we have taken uik = 0.495 GeV and 
fx = 0.1598 GeV. In the chiral limit (mp = 0) we use the "experimental" 

9 Since it is known that the scale is poorly fixed from the nucleon mass, we have 
not taken this estimate into account. 
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value f x = 0.1282 GeV obtained by considering fx and f n , as functions of m 2 K 
and m% respectively(PCAC) and by extrapolating them naively to vanishing 
m P . 



5.2 Standard lattice measurements of the dimensionless chiral condensate 



We will now present several "standard" measurements of the subtracted chiral 
condensate based on eqs. (29), (30) and (33). In theory all three equations 
are equivalent in the continuum limit. However, since they are implemented 
in practice at fixed lattice spacing (fixed (3) and non-zero quark mass, they 
are subject to systematic errors. Thus, by obtaining the chiral condensate in 
several ways, we are able to study the systematic errors introduced by the 
finiteness of the lattice spacing, the extrapolations to the chiral limit etc. 

Once the hadronic correlation functions have been fitted and the lattice masses 
and matrix elements have been extracted, we must extrapolate the results to 
the chiral limit. The standard way consists in determining the chiral point kc 
by fitting the pseudoscalar mass Mj> as a linear function of the quark masses 

Mp = C HS (-- — ) . (36) 

\K K C J 



The slope C HS related to the hadron spectrum (HS) is also determined from 
this fit. The problem is that this relation is based on eq. (35) which, for 
our Clover dataset is only satisfied to 0(a) and therefore it should not be 
implemented for this action (see Appendix C). In order to avoid this problem, 
we fit instead ap as a linear function of Mj> 

2a P = 7^wl M P > ( 37 ) 



and determine the slope C AWI related to the Axial Wis (AW I). We have 
verified for the Wilson data that ap vanishes at k c within our statistical 
errors. 

The most straightforward way of computing the subtracted condensate is that 
of eq. (29), where the term (m — m) has been computed from ap, with the 
aid of eq. (25): 

= J™ ap^(V(x)V\0)) . (38) 

Z A M^O V 



Unfortunately this method is unreliable: at (3 — 6.0 and 6.2 we obtained 
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results which are clearly incompatible with those obtained from other methods; 
at /3 — 6.4 the statistical error overwhelms the signal. The reason is that 
the correlation function (V(x)V^(0)), integrated over all space-time, has a 
contact term at x xs which from dimensional counting is seen to behave as 
I /a 2 . Although its contribution to the chiral condensate vanishes in the chiral 
limit (it is to be multiplied by (mo — fn)), it is the dominant contribution 
at the values of lattice spacing and quark mass of the simulation. Therefore 
the extrapolation to m c becomes unreliable and with big statistical errors. 
Thus, this computation of the condensate must be discarded. This effect was 
first observed (with much poorer statistics and control of the systematics) in 
ref. [12]. However we note, that in ref. [13] a variant of the above determination, 
performed with smeared operators, gave satisfactory results in the chiral limit. 
Presumably, smearing reduces the influence of the contact term on the data. 

Two robust ways of computing the subtracted chiral condensate, which avoid 
contact terms, are based on eqs. (30) and (33). Once more eq. (25) is used to 
compute the term (m — m). We find for (xx) : 



(xxh 

Z A 

(xxh 

Z A 



17777 lim 

2C AWI Ml^o 



(0\V(0)\P) 



(39) 



= — C 



AWI 



lim Ft 



An- 



other computations of the subtracted chiral condensate are also possible by 
using eqs. (24), (35) and (36) in order to express (m —m) in terms of (mo-mc). 
We avoid this option for two reasons. First, it requires the ratio Zp/Z$ which, 
compared to Za, has a larger error (see results in Table 2). Second, it involves 
eq. (35) which cannot be implemented with our Clover data. 



Run 


(xx)i/Z A 


(xxh/z A 


a 3 (xx)i 


a~ 3 (xx)2 


C60a 
C60b 
W60 


0.0040(3) 
0.0039(2) 
0.0045(3) 


0.0040(3) 
0.0039(2) 
0.0051(4) 


0.044(4) 
0.045(3) 
0.048(4) 


0.044(4) 
0.046(3) 
0.053(6) 


C62 
W62 


0.0018(2) 
0.0019(1) 


0.0015(2) 
0.0016(1) 


0.041(7) 
0.046(5) 


0.034(5) 
0.038(4) 


C64 
W64 


0.0007(1) 
0.0008(1) 


0.0007(1) 
0.0007(1) 


0.047(12) 
0.053(11) 


0.049(8) 
0.047(7) 



Table 3 



The subtracted chiral condensate determined from eqs. (39) in lattice units (second 
and third columns) and in physical, GeV 3 , units (fourth and fifth columns). An 
overall minus sign has been omitted from the data. All errors are statistical. 
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In Table 3 we present our results. Our comments are the following: 

• The results from runs C60a and C60b, which only differ by a volume factor 
of about 1.5, are perfectly compatible. We see no finite volume effects. 

• All (xx)i/Za and (xx)2/Za results, obtained with the Wilson action at the 
same (3, are compatible within two standard deviations. Those from the 
Clover action are compatible within one standard deviation. This suggests 
that our chiral extrapolations are robust, especially in the Clover case where 
the discretization error is reduced. 

• The (3 = 6.4 results have a larger relative error for both actions. This is 
probably due to the small physical volume at this value of (3. We recall that 
the quark mass measured at this coupling on the same dataset in ref. [9] 
also suffered from finite lattice effects. Thus the (3 = 6.4 results will not be 
taken into account. 

Earlier works also computed the subtracted condensate at (3 = 6.0 with the 
Wilson action. The methods used were all based on eq.(33) (GMOR), so they 
are comparable to our {xX)2/Za value: 



We see that our result is compatible with previous determinations, but the 
statistical accuracy is greatly improved. 

The standard method for computing the chiral condensate would now involve 
renormalizing (xx)i or (XX)2 by multiplying it with Zp (c.f. eqs.(27) and (28)) 
and expressing it in physical units by multiplying it with a~ 3 (c.f. eq.(34)). The 
latter step is the source of error amplification, due to the cubic power of the 
lattice spacing. Our 8% overall estimate of the error due to the lattice calibra- 
tion would be amplified to 24%. This increased error is a serious shortcoming 
of all standard methods of computing the condensate on the lattice. 



5.3 Lattice measurements of the physical chiral condensate 

We will now present a better method of computing the chiral condensate on 
the lattice. Rather than the standard two-step approach (i.e. first measure 
the dimensionless condensate in terms of dimensionless quantities and then 
multiply it by a~ 3 ) we will obtain it in physical units, by using the physical 
value of the pseudoscalar decay constant multiplied by a -1 . Thus we write the 



(XX)2/Z A = -0.0051(4) 
(xx)/Za = -0.006(1) 
(XX)/Z A = -0.006(3) 



Ref. [12] 
Ref. [13] . 



This work 



(40) 
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GMOR relation (33) for the renormalized condensate as follows 

= -\a' l f x Z s C HS . (41) 

Use of eqs.(24) and (36) has been made. Note that f x = 0.1282 GeV is the 
"experimental" value in physical units (see subsect. 5.1). The computation 
of the condensate based on the above formula has several advantages. The 
most important is that, by expressing it in terms of f x , we are left with only 
one power of the UV cutoff a -1 . Thus, we avoid the error amplification of 
the standard methods. Another advantage is that there is no extrapolation 
to the chiral limit. Instead, we only need to determine the slope C HS of the 
squared mass of the pseudoscalar meson with respect to the quark mass. This 
is reminiscent of similar procedures used in spectroscopic studies in refs. [8,33]. 
Note that since we work at k values typical of the strange quark mass, we are 
implicitly assuming that the slope will not change in the chiral region. There 
is no way round this assumption, besides simulating at lighter quark masses. 
Finally, another advantage of eq. (41) is that the renormalization constant is 
now Zg, which is much better determined than Zp with the NP method. 

The only drawback of the above determination is that it cannot be used on our 
Clover data as explained in Appendix C. The alternative is to write eq. (33), 
in terms of eq.(37), as follows 

J-(^) 2 = -\ a -ifl?Lc AWI . (42) 
N f 2 Jx Z A y ' 

Also this determination has the advantages of having only a single power of 
the inverse lattice spacing and no explicit chiral extrapolation; a disadvantage 
is the larger systematic error of Zp. 

Since we use the NP estimates of the renormalization constants in both de- 
terminations of eqs. (41) and (42), our results are obtained in the RI scheme 
at a renormalization scale a -1 . Subsequently, perturbation theory is used to 
NNLO to express them in the MS scheme (c.f. eq. (12)) and run them to 
the conventional reference scale of 2 GeV (c.f. eq. (13)). In so doing, two 
further uncertainties are introduced. The first is the 0(a^) error due to the 
truncation of the NNLO perturbative series which we consider negligible. The 
second uncertainty arises because the initial scale \i ~ a -1 is known with an 
8% precision. Since the dependence of the perturbative coefficients on the scale 
is logarithmic, we find that the error introduced is at most 2% and will be 
ignored. 

In conclusion there are three important errors in our results. Firstly, we have 
a purely statistical error of the lattice quantities C HS and C AWI . Secondly, 
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there is the 8% error of statistical and systematic origin from the lattice cali- 
bration (the a~ l factor in eqs.(41) and (42)). Finally, we have the error in the 
determination of the renormalization constants (c.f. Table 2) which is also of 
both statistical and systematic origin. 

In table 4 we present our results for the chiral condensate in the MS scheme 
at 2 GeV. Our comments are as follows: 

• At fixed (3, the two estimates from (ipip)^ 15 and (iptp)^ obtained with the 
Wilson action are perfectly compatible within the statistical errors. 

• At (5 — 6.0, the {iJji/j)^ s results obtained with the Clover and Wilson actions 
are perfectly compatible within two standard deviations of the statistical 
error; at /3 — 6.2 they are compatible within one standard deviation of the 
statistical error. This implies that discretization effects are under control. 

• Even within our increased statistical accuracy it is not possible to see a 
systematic dependence of the chiral condensate in the f3 range 6.0-6.2. The 
results at /3 — 6.4 apparently show such a dependence, but, as we have 
argued in subsect. 5.2, this is more likely to be a finite size effect. Thus, we 
will again not take into consideration results at /3 — 6.4. 

• As previously stated, the systematic error due to the NP renormalization 
constants is small in all {$vjj)Y s results and dominant in the {ip7p)^ s results. 
However we cannot exclude that the very small error in Z$ is underesti- 
mated. A conservative attitude consists is regarding the larger error of Zp 
in the {i>ip)2 LS results as more realistic. 

In order to appreciate the effects of the NP renormalization on the chiral 
condensate, we now present the results at /3 — 6.2 obtained from the per- 
turbative values of Z$ and Zp. We use the results of Boosted PT (BPT) 
(see refs. [19,28]). For the Clover action we obtain (V^)F^ = -0.0185(10), 
whereas from the Wilson action we have (•0*0)f f5 = —0.0158(5) and ( , ipip)^ IS = 
—0.0205(7). It is obvious that the excellent compatibility of the corresponding 
results, obtained with NP renormalization constants, has been destroyed by 
PT. 



5.4 Final result 

From the discussion in the previous section on the results of Table 4 we can now 
derive our best estimate for the chiral condensate. Our best results are those 
at (3 — 6.2 (high enough UV cutoff without finite volume effects). Although 
all results at this (3 value are compatible, the Clover estimate is in principle 
preferable, as it has O^a) discretization errors. Thus we give as our best 
estimate the C62 result of Table 4: 
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Run 


C HS 




qAWI 




C60a 
C60b 
W60 


2.98(8) 
3.04(7) 
2.40(5) 


0.0150(3)(2) 


3.9(1) 
4.1(1) 
3.01(7) 


0.0141(5)(21) 
0.0146(3)(21) 
0.0152(4)(20) 


C62 
W62 


2.9(1) 
2.52(8) 


0.0156(5)(2) 


3.7(2) 
2.98(9) 


0.0147(8)(16) 
0.0158(5)(16) 


C64 
W64 


3.5(2) 
2.9(1) 


0.0172(7)(2) 


4.2(2) 
3.2(1) 


0.0187(9)(10) 
0.0179(8)(13) 


Table < 


1 



The pure lattice quantities C HS , C AWI and the corresponding chiral condensates 
(in GeV 3 ) at a scale fi = 2 GeV, in the MS scheme. The first error is statistical, the 
second comes from the NP renormalization constants. An overall minus sign has 
been suppressed in the results of the chiral condensate. 

-^r(^) m (fi = 2GeV) = - (0.0147 ± 0.0008 ± 0.0016 ± 0.0012) GeV 3 
N f 

= - [(245 ± 4 ± 9 ± 7 MeV)] 3 , (43) 

where the first error is statistical, the second is due to the NP renormalization 
and the third due to the lattice calibration. For the first time we are able to 
estimate the systematic uncertainties in a reliable way. We have also seen that 
our analysis is in accordance (with improved statistical accuracy) with earlier 
lattice results [12,13]; c.f. eq. (40). The results cited here disagree by a factor 
of 2 with the recent lattice computation of ref. [34]. Similar discrepancies were 
seen between the quark mass estimates of this work and, say, ref. [9]. The 
two effects are not independent, since the GMOR relation connects the chiral 
condensate to the quark mass. Presumably, the results of ref. [34] are flawed 
by unreliable extrapolations of the lattice data obtained at several UV cutoff 
values to the continuum limit. 

Our result is also in agreement with the one obtained from sum rules (see [14] 
and references therein): 

-^{W>) m (A* = 2 GeV) = - (0.014 ± 0.002) GeV 3 
N f 

= - [(242 ± 9) MeV] 3 . (44) 



For completeness we also give the quenched RGI result denned in eq. (14): 
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_L(^) RGI = - (0.0088 ± 0.0005 ± 0.001 ± 0.0007) GeV 3 
N f 

= - [(206 ± 4 ± 8 ± 5 MeV)] 3 . (45) 
The errors correspond to those of eq. (43). 

The reference scale \x = 2GeV is the standard lattice choice. In order to facil- 
itate comparison with other determinations of the condensate, conventionally 
given at the scale \x = lGeV, we run our results to the latter scale, using 
NNLO perturbation theory. We find 

-^ r (4>'ip) m (// = lGeV) = - (0.0124 ± 0.0007 ± 0.0014 ± 0.0010) GeV 3 
N f 

= - [(231 ± 4 ± 8 ± 6 MeV)] 3 , (46) 
to be compared to the sum rule result of ref. [14]: 

-^r{W) m (/i = lGeV) = - (0.012 ± 0.002) GeV 3 
N f 

= - [(229 ± 9) MeV] 3 . (47) 



6 Conclusions 

We have computed the QCD chiral condensate from first principles in the 
framework of the lattice regularization with Wilson fermions. Our result has 
been obtained in the quenched approximation. We have been particularly care- 
ful in understanding and controlling all important sources of error (except for 
quenching): (i) large field configuration ensembles have been generated in or- 
der to minimize the statistical error; (ii) two actions (Wilson and tree-level 
Clover) and several gauge couplings have been used to control the discretiza- 
tion error; (iii) two lattice volumes (at one gauge coupling) ensure some con- 
trol of the finite size effects; (iv) several methods of extracting the chiral 
condensate, equivalent only in the chiral limit, give us some confidence about 
the extrapolations to the chiral limit; (v) all necessary renormalizations have 
been performed non-perturbatively to avoid large tadpole contributions which 
distort lattice perturbative renormalization at 1-loop; (vi) the RI - MS renor- 
malization scheme matching and the RG running from the lattice scale a -1 to 
the usual scale 2 GeV have been done in continuum perturbation theory at 
NNLO. Our final estimate for the chiral condensate is given in eq.(43). It is 
compatible with previous lattice determinations with a much smaller statisti- 
cal error and a careful estimate of the systematic error. It is also compatible 
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with the sum rule result. Thus we conclude that the chiral condensate is an 
order parameter of the spontaneous breaking of chiral symmetry in QCD with 
massless quarks. 
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APPENDICES 



A The RI renormalization scheme 



We shortly review the NP method for the renormalization constants of quark 
bilinear operators. The full discussion of the method can be found in ref. [25]. 
The most up to date results, which we use in the present work, can be found 
in ref. [28]. 

Given a quark bilinear 0[ = ipT(\^ /2)ip (T is any Dirac matrix), we consider 
the operator insertion in the quark propagator 

Go(p) = /d^expfefe - ^MsOOfcO)**,)) , (A.1) 

from which we obtain the amputated Green function computed between off- 
shell quark states of momentum p in the Landau gauge: 

Ao(pa) = S(paY l G (pa)S(pa)- 1 . (A.2) 



The above quantity is computed non-perturbatively via Monte Carlo simu- 
lations in the Landau gauge. [25]. Any effects from Gribov copies are small 
[24]; spurious solutions [35] have not been considered. The renormalization 
Z o (na,g ) of Or is determined by the RI condition 



Z (i^a)Z 1 (//a)Tr P Ao(pa) 



(A.3) 



where Po is a projector chosen so that the above condition is satisfied at tree 
level [25] and Z q is the wave function renormalization. Since the projected 
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amputated Green function can be computed non-perturbatively with Monte 
Carlo simulations, eq.(A.3) can be solved for Z Q at a fixed scale \i and cutoff 
a provided that Z q is known. 

In order for the RI scheme to be compatible (at large scales //) with the Wis, 
Z q must be defined as [25] 

(A.4) 

p 2 =H 2 



The above definition is not convenient, as it requires differentiation with re- 
spect to the discrete variable p. Thus, following ref. [25], we have opted for 
the definition 



Z' q (fw) = - 



. 1 TrT,ij,=i,A^nSm(p^a)S(pay 



12 4 4 sin (p^a) 



(A.5) 



which, in the Landau Gauge, differs from Z q by a finite term of order a 2 . The 
matching coefficient can be computed using continuum perturbation theory. 
From refs. [28,30] we quote up to order a 2 , 



Z' q (47T) 



* = 1--^A® + ... , (A.6) 



where 



Af = ^L— l l (3 + 22N 2 - AN c N f ) (A.7) 



16 N 2 



(Nf = in the quenched approximation). Note that at a typical scale fi ~ 
2 GeV, A q represents a tiny correction, smaller than 1%. The renormalization 
constants used in the present work have been corrected by this factor (c.f. 
[28]). 



B RI/MS matching coefficient and RG evolution coefficient 



The matching coefficient AZ RI / MS which transforms the renormalized chiral 
condensate from the RI to the MS scheme is that of the renormalization con- 
stant Z P (or Zs). It has been calculated in PT to 0(a 2 s ) (see refs. [28,30] for 



26 



details). The result is 



az ri/ms = x + ?M c {i) + ?^ C ^ + 0{al) , (B.l) 
An (An) 2 



where 



C^^l (B.2) 



n_ 8(iv c 2 -i) 

AN C 

= ~^ (-309 + 3029 N 2 C 
yo iv c 

-288 C 3 - 576 N 2 C ( 3 - 356 iV c N f ) 

and Cs = 1.20206 



The evolution coefficient of the chiral condensate is derived from a standard 
RG analysis. Its running is determined by the evolution of the strong coupling 
a s (/i) (i.e. the Callan-Symanzik (5 function) and the scalar operator (i.e. its 
anomalous dimension 75 (a s ). In any scheme which respects vector symmetry 
(i.e. the vector Wis are valid) the renormalization of the scalar operator is the 
inverse of the renormalization of the quark mass 10 (see for example eq.(17)). 
Thus, we can readily use the results of refs. [36] obtained in the MS scheme 
(with NDR regularization) at the NNLO. The Callan -Symanzik (3 function 
and the quark mass anomalous dimension are given by 



An dfi 2 \An) ' 

^•)=- 2 ^ 2 $ = g^(S)' +1 . (b.3) 



where 



/3o=yiV c -| J V / (B.4) 

fl 3 c 3 c f N c i 

/J p = 2857 W-lf N _ 205 , _ , 

54 c AN? i 36 V c I s 

70 

r2 



U1 Kr2, T , H(iV c 2 -l) Ar2 , 79 



N'Nf + — ±-2 -N f + — N c N f 

54 c / 18 N r f 54 f 



10 



More generally, the anomalous dimensions of the two-fermion operators of eqs. (5) 



and (6) satisfy jv = 1A = and 7s = 7p = — 7 m - 
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Im 



Im 



^(2) 
Im 



N 2 - ' 

N c 

N 2 



'23 

+N f [ -N c 



3 203 
~4 + ~Y2 



129 



N 2 C - jjiV^) 



129 Ar2 11413 Ar4 



1177 *r3 

54 c 



108 



35 



(B.5) 



Note that the scheme dependence settles in only at the two-loop order for 
7 m (a s ) and three-loop order for f3(a s ). Finally, to NNLO the evolution coeffi- 
cient of the chiral condensate is [36] 



B. 



1 / a s (/i) 

2 I 4tt 



with ft = andjfc = 7sV( 2 A>) = -7^/(W- The running coupling 

Oi s (n) at NNLO in MS scheme is given by 



«f% 2 ._ 1 A In ln(g 2 ) 

4tt W j /5 ln(g 2 ) (31 hiV) 

+ ^]^y (A 2 In 2 ln(g 2 ) - (3f In ln(g 2 ) + /3ffo - /3 2 ) , (B.7) 

where q 2 = (fx/ Aqq D ) 2 . Since we work in the quenched approximation Nf = 
in all formulae. The QCD scale in the quenched approximation has been set 
to Agg D = 0.251 ± 0.021 GeV; see ref. [37]. 



C Lattice Wis with Clover fermions 



In this Appendix we discuss the modifications which are necessary for a compu- 
tation of the chiral condensate with the tree level Clover improved action. This 
involves the discussion of several technical details. At tree-level, the Clover ac- 
tion is defined as follows [10]: 

S c = S f - a 4 ^2 ji9o^(x)a IM/ F lu/ (x)ip(x) , (CI) 
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with F lll/ (x) the clover- leaf discretization of the field tensor. At this order all 
Oiag^ 1 ln n a) terms, which are effectively of 0(a) in the scaling limit (g$ ~ 
1/lna), are eliminated from correlation functions. At leading- log level the 
improvement of local operators can be expressed as a rotation of the fermion 
fields [11]: 

0*.(x) = i) R {x)V^ R {x) , (C.2) 



where the rotated fields are defined through 



1 - 



a p 



if>(x) 



ip R (x) = ip(x) 



a Jp 



(C.3) 



and p and p are symmetric lattice discretizations of the covariant derivatives 
(see ref. [11] for their definition). The improved operators O 1 differ from the 
original ones by terms proportional to the cutoff. When O 1 is inserted in a 
Green function, these extra terms combine with the a -1 UV divergences to 
give finite contributions. Consequently, Clover renormalization constants differ 
from the Wilson ones by finite terms. 

We denote by S^x — y) — (ip(x)ijj(y)) the quark propagator obtained by 
solving the Dirac equation of the Clover action. The tree- level 0(a) improved 
quark propagator is {■ifj R (x) / ijj R (y)); in terms of S 1 (x — y) it is given by (see 
ref. [19] for details): 



(1> R (x)$ R (y))=S eff (x 



y) + -S(x-y) + 0(a 2 



(C.4) 



where the effective rotated propagator S e ^(x — y) is defined as 



S eff {x-y) = 



1 - g P (*) 



S\x-y) 



1 + 2^0/) 



(C.5) 



This is the Clover propagator we use in our computations. The reason it 
contains an 0(a 2 ) term in its definition is that it is readily computable in 
numerical simulations [38] . Thus, it has been extensively used in several Clover 
improved lattice QCD computations of on-shell matrix elements. In computing 
correlation functions between external on-shell hadron states, the 5-function 
never contributes. To illustrate the point, consider a two-point correlation 
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function of two bilinear improved operators Oj(x) = ip R (x)rjip R (x) (where 
j = 1,2) given by 



(0[(x)Oi(y)) = {0[{x)Oi{y)) n4ct (C.6) 
~5{x - y)Tr [r!r 2 5 e// (a; - y)} - ^5(x - y)Tr [T 2 T 1 S eff (x - y) 



where 



{0[(x)0{(y)) n4ct = -(Tr [s°"(x - y)V 2 S^{y - x)^]) (C.7) 



(i.e. the subscript nlpct stands for "no Jf) contact terms"). Note that the nlpct 
correlation functions are traces of the effective propagator S e ^(x — y). How- 
ever, in the case of the chiral condensate, c.f. eq. (29), the space-time points 
of the two-point correlation function are allowed to coincide. The on-shell ma- 
trix elements argument is then invalidated and the complete contribution of 
eq. (C.6) must be used with the 5-function giving rise to contact terms. 

The above improvement scheme is by no means unique. Other choices of quark 
field rotations are also possible, as discussed in refs. [19,22,39]. The choice 
made here is not necessarily the most convenient for improving lattice Wis. 
As previously pointed out, it is however, the one implemented in all APE 
simulations with tree-level Clover improvement. Thus, we need to examine 
how (9(a)-improved Wis can be obtained for this choice of rotation (i.e. this 
choice of improved operators). 

Upon performing an axial variation of the improved pseudoscalar operator 
P l9 (0) with the theory defined in terms of the Clover action, we obtain 



s(P l9 (Q)) , t 

Sa^ A (x) 

*<f^f > = fl4 E^{4(x)P ;9 (0))-a 4 2 mo (P / W^(0)) (C.8) 



-a\Xf{x)P l9 ({])) 



This is eq. (18) with the difference that the variation of the Clover term in the 
action modifies the expression for X* by a term proportional to a^a^F^ip 
(c.f. eq. (C.l); F M „ is the dual field tensor). The above WI is not 0(a) improved 
for three reasons: (i) the correlation functions contain unimproved operators; 
(ii) the asymmetric derivative is not improved (V£ = d% + 0(a)); (iii) the 
variation of the improved operator P l9 (the l.h.s. of the above) does not yield 
the vacuum expectation value of the improved scalar density (S IO (0)) (c.f. 
eq. (18)). The conclusion is that improved Wis cannot be obtained by the 
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standard procedure of performing symmetry transformations on improved op- 
erators and the Clover improved action. Instead, we simply impose the validity 
of Wis (such as that of eq. (18)) with symmetric derivatives and correlation 
functions of improved operators. Since all equations of sect. 4, implemented at 
fixed inverse coupling /3, suffer from 0(a) corrections, we can add by hand the 
0(a) terms which render the correlation functions improved. Consistency re- 
quires that the mass subtraction m(mo) is also modified by 0(a) terms so that 
(m — m / )P / ^ is O(a) improved. Consequently, the mass difference (m — m 1 ) 
vanishes at a critical point tUq. In other words, we define 2// as: 



where V^, is the symmetric lattice derivative. This equation, with Z A i com- 
puted, say, from improved Wis [19,24,40] is a definition of the improved quan- 
tity (m — m 1 ). With this definition, this mass difference is renormalized by 
Z~j (the improved version of eq. (23)). 

We can now write down the improved WI for the chiral condensate. It is the 
WI (29) expressed in terms of improved operators and masses. Also, contact 
terms must be properly taken into account according to eq. (C.6). We obtain 



Apart from these modifications, all other expressions for the chiral condensate 
(see sect. 5) remain valid, provided all quantities are improved. 

Similar arguments can be used in the case of vector Wis. It is not possible 
to obtain Clover-improved Wis by the standard vector variations of improved 
operators. Instead, we impose the validity of WI (16) with all correlation 
functions expressed in terms of improved operators. This implies an improved 
definition of the subtracted mass (mo — uiq), so that the renormalization con- 
stants Z s i and obey eq. (17). Moreover, as stressed repeatedly in this 
work, the simple relationship between the quark mass and the hopping pa- 
rameter of eq. (35), satisfied to 0(a), needs to be improved. As this is not a 
straightforward task with the rotated fields of eq.(C.3), we have avoided use of 
eq.(35) when deriving our Clover results. This point has been explained also 
in ref. [9], where an improved version of eq. (35) valid for a different rotation 
of the quark fields is given. 



2/ = 



2(mo — m 



:') _Jd 3 xV° x (A I f (x)P I f(0)) 



(C.9) 



Z A i 



Jd 3 x(P I f(x)P I f(0)) 
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A final word of caution is in place here. The relationships (17) and (23) con- 
strain the renormalization constants involved, so that the Wis are recovered in 
the continuum. These constraints relate not only the anomalous dimensions, 
but also the finite parts of these renormalization constants. The determina- 
tion of the finite part of a renormalization constant depends on the choice of 
lattice operator^ and renormalization condition. There are choices for which 
the Wis are not satisfied, and consequently eqs. (17) and (23) are violated 
by finite terms. For example, in 1-loop perturbation theory with the Clover 
action, has been obtained by renormalizing the Clover quark propagator 
(ipip) (in the MS scheme) [21]. Compared to the result for Z s i obtained in 
ref. [22], we see that eq. (17) is violated by the finite term arising from the 
quark field rotations in S 1 [9]. Agreement would have been reached, had the ro- 
tated propagator (ip R ijj R ) been used to obtain Z^. The discrepancy is not due 
to the choice of scheme, but to the choice of an unrotated quark propagator. 
On the other hand, the renormalization constants used in this work, obtained 
in the RI scheme with the NP method of [25], are perfectly compatible with 
the Clover Wis. 
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